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A MATHEMATICAL EXAMINATION OF THE PRESS MODEL 


FOR ATMOSPHERIC TURBULENCE 

Kenneth Sldwell* 

Langley Research Center 

SUMMARY 

The random process used to model atmospheric turbulence in aircraft response 
problems is examined on a mathematical basis. The first, second, and higher order 
probability density and characteristic functions are developed. The concepts of the 
Press model lead to an approximate procedure for the analysis of the response of lin- 
ear dynamic systems to a class of non-Gaussian random processes. The Press model 
accounts for both the Gaussian and non-Gaussian forms of measured turbulence data. 

The nonstationary aspects of measured data are explicitly described by the transition 
properties of the random process. The effects of the distribution of the intensity pro- 
cess upon calculated exceedances are examined. It is concluded that the Press model 
with a Gaussian intensity distribution gives a conservative prediction of limit load values. 

INTRODUCTION 

The random process developed by Press and his associates at NASA (refs. 1 to 3) 
is used to model the random nature of atmospheric turbulence in many aeronautical appli- 
cations. The process is an application of random process theory to the problems of air- 
craft response to atmospheric turbulence. It was developed to account for the exceed- 
ance properties of atmospheric tubulence indicated by data from VGH recorders on many 
aircraft. The turbulence was considered to consist of a series of local regions or 
"patches," each having a stationary Gaussian distribution of turbulence velocity compo- 
nents with a fixed intensity or standard deviation. The intensity was then varied by a 
separate random process in order to account for the apparent nonstationary and non- 
Gaussian structure of the atmospheric turbulence. A schematic example of the pattern 
of locally Gaussian regions with intensity variations between the regions is shown in 
figure 1, 

The probabilistic structure of the resulting random process was first examined by 
Bullen. (See refs. 4 and 5.) Other studies of the random process are the reports by 
Dutton and Thompson (ref. 6) and by Reeves (ref, 7). The application of the process to 

*This research was accomplished while the author held a National Research Coim- 
cil Postdoctoral Resident Research Associateship at NASA Langley Research Center. 



atmospheric turbulence data has also been examined by Coupry (refs. 8 and 9) who refer- 
red to the development as the Press process, a terminology which is followed in the 
present report. The associated exponential type of exceedance expression for the aver- 
age rate of exceeding a level of aircraft response has been the basis for the evaluation of 
an extensive amount of experimental data. Some of this work has been discussed in ref- 
erences 10 to 13. Several studies have developed the associated analytical techniques 
(the "power spectral density method") for application in aircraft strength design. (See 
refs. 14 to 16.) These procedures are included in some versions of aircraft strength 
criteria. (See refs. 17 and 18.) 

The present study examines the mathematical basis of the Press process and exam- 
ines several applications both to the measurement of atmospheric turbulence and to the 
associated problems of aircraft response. The primary conclusion is that the basic con- 
cepts of the Press process, which were originally developed largely on an intuitive basis, 
can be expressed with an explicit mathematical model. The study is divided into two 
parts. The first part is a general mathematical formulation of the Press process. Ini- 
tially, the first order distributions of the process are developed, that is, distributions of 
an ensemble at a single time point. Then the second and higher order distributions are 
developed; this development leads to a discussion of the dynamics of the random process. 

The second part of the study considers several applications of the Press process. 
Two specific mathematical problems are discussed. First, the validity of currently used 
anal 3 dical procedures as a method for the analysis of dynamic system response to a class 
of non-Gaussian random processes is discussed. Second, the ability of the Press process 
to account for the nonstationary behavior of atmospheric turbulence is discussed. It is 
also shown that the transition properties of the Press process explicitly account for both 
the Gaussian and non-Gaussian forms of experimental data. The use of the Press pro- 
cess as a specific model for atmospheric turbulence in aeronautical applications is then 
discussed. The effects of the choice of the intensity process upon both measured and 
calculated exceedance data are considered, 

SYMBOLS 

A standard deviation of R process 

b standard deviation of S process 

bj intensity parameter 

C( ) characteristic function of subscripted process 
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det 

{dR( )} 
E[ ] 

F( ) 

H 

Hij 

h(r) 

i 

inv 

KnO 

M 4 


determinant of a matrix 

spectral process associated with R process 
ensemble average 

spectral distribution function of subscripted process 
covariance matrix of vector R process 
elements of covariance matrix of vector R process 

impulsive response function of a linear system 
unit of imaginaries, \f^ 
inverse of a matrix 

modified Bessel function of order n (ref. 22) 

flatness factor (ratio of fourth moment to square of second moment) 


m 

N( ) 
No 


Pi 

P() 

Pc(z|s) 


dimension of vector process 


expected rate of positive slope crossings of indicated level 


1 

expected rate of positive slope zero crossings, 

Z7T 


index of generalized Press process (=1/2, 1, 3/2, . . .) 


probability parameter 


probability density function of subscripted process 

conditional probability density function of z process conditional on value s 


r,R random process 
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I 


s,S random process ("intensity") 

T time interval 

t time 

Ucr turbulence or gust intensity parameter 

x, X quasi -steady random process (Press or generalized Press) 

xq mean value of x 

y, Y generalized Press random process 

z, Z product or Press random process 

a = \/l - p2 

/3 Ab0 

r( ) gamma function (ref. 22) 

6ij Kronecker 6-S3mibol 

9 Fourier transform variable for characteristic functions 

p correlation coefficient, correlation function, or autocorrelation function 

S,a random process 

T time -difference variable 

#(o)) power spectral density function of subscripted process, -« g w s oo 

^( ) autocovariance function of subscripted process 

CO frequency, Fourier transform variable of t 
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time average 


<> 

Superscript: 

a integer index 

Subscripts: 

c conditional 

d derivative 

i integer index 

in input random process 

j integer index 

I limit 

out output random process 

tr transition 

A bar over a S5mibol denotes a vector. Dots over symbols denote derivatives with 
respect to time. Primes denote derivatives with respect to the argument. A tilde over 
a symbol denotes a variable of integration. 

FIRST ORDER DISTRIBUTIONS 

The formulation of the Press process is developed in the present section. The 
development extends some previous work which has been done on the mathematical form- 
ulation of the Press process in references 4 and 6. The basis of the present development 
is the interpretation of the Press process as the product of two independent random pro- 
cesses, a concept first stated by Reeves (ref. 7). 

The discussion is restricted to first order distributions; that is, the process is 
considered at a single time point only. Thus, the treatment of the time evolution of the 
random process is not considered in this section. The development follows standard 
analytical techniques of probability theory; appropriate references are the texts by 
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Papoulis (ref, 19) and by Lin (ref. 20). All the Indicated statistical moments are assumed 
to exist. All the random processes are specified to have zero mean values. 

Product Process 

The random process is defined by the product of two component processes, 

z(ti) = r(ti) s(ti) (1) 

The two component processes, R and S, are specified to be statistically independent. 

The probability functions of the product process Z are completely determined by 
the functions of the two component processes. The joint probability density function of 
the product process and one of its components is determined by the product relation 
(eq, (1)) and the independence of the two component processes 

Pzs(z,s) = Pp(r = z/s)Ps(s) (2) 

The probability density function of the product process itself is obtained from the joint 
density function by using the consistency relation, which determines the marginal distri- 
bution from the joint probability density function 

Pz^^) 7^ Pr(^ = z/s)pg(s) ds (3) 

'^-oo I® I 

The R component process can be related to an associated conditional probability 
density function. By using the definition of conditional probability and the form of the 
joint probability density function (eq. (2)), 

Pc(z:|s) = Pr(r = z/s) (4) 

The probability density function of the product process can be written in terms of the 
conditional probability density function as 



The probability density function of the product process can be develc 5 )ed from 
either equation (3) or equation (5), depending upon whether the R component process is 
specified either directly or indirectly by the associated conditional process. The his- 
torical approach (ref. 1) used the idea of the conditional probability and developed equa- 
tion (4) through the specific form assumed for the conditional process. The interpreta- 
tion of the resulting random process as a product process was introduced by Reeves 
(ref. 7). The product relation (eq. (1)) is particularly useful in the subsequent develop- 
ment since it leads directly to some of the important properties of the Press process and 
its extensions. 

The product process can also be developed in terms of its characteristic function 


Cz 



e^®^ Pz(z) dz 


( 6 ) 


The characteristic fxmction is developed from equation (3) by using the properties of the 
associated Fourier transformations (ref. 21) 


Cz 



Cr(0s) Pg(s) ds 


( 7 ) 


The formulation of the product process in terms of characteristic functions is conveni- 
ent since it usually leads to simpler functional relations. 

The statistical moments of the product process can be developed from those of the 
component processes; either from both the probability density and the characteristic 
functions or directly from the product relation and the independence of the component 
processes 


E[z«(ti)] = Ep(ti)]Ep(ti)] (8) 

The preceding development considers the product process in general terms. The 
subsequent development is restricted to the specific forms of the two component pro- 
cesses which are useful in applications as a model for atmospheric turbulence. 

Press Process 

The original development (ref. 1) of the Press random process specified the R 
process to be Gaussian under the locally Gaussian condition. All the present develop- 
ment is restricted to this case. The reason for this choice is discussed later in the 
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section on the quasi -steady approximation. Several distributions were considered for 
the S process in reference 1. A Gaussian distribution was selected, based on the cor- 
respondence between the resulting exponential exceedance expression and measured 
atmospheric turbulence data. Almost all the subsequent applications of the Press model 
have been based upon the S component being a Gaussian random process. This partic- 
ular case is referred to as the Press process. 

Scalar-scalar. - The Press process in scalar form is generated by the product of 
two scalar Gaussian independent component processes, R and S. 


[r2(ti)^ 

= a2" 

> 

s2(ti) 

= b2 

J 


The statistical moments of the Press process can be determined directly from the 
product relation and the independence of the R and S process. 


E 



= 0 


E 


z^(ti) = [l 


3 • 5 


. m - 1 )] 


2A^b^ 


(e, odd) 

!> ( 10 ) 

(f, even) 


The probability density functions of the component processes have the Gaussian functional 
form: 


Pr(r) 



Ps(s) 



(-00 = r = oo) 




(-0O = S = oo) 


( 11 ) 




This notation is fairly standard in the aeronautical literature, except that the S compo- 
nent process is usually restricted to nonnegative values. However, the S process must 
be defined with both positive and negative values in order to avoid some mathematical 
problems in the subsequent formulation. This approach does not lead to any incompati- 
bility with the traditional approach to the Press model. The possible restriction to non- 
negative values is based upon the historical development of the Press process through 
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the concept of a conditional process, using equations (4) and (11), where the conditional 
standard deviation is the product of the standard deviation of the R process and the 
modulus of the value of the S process. In order to satisfy this interpretation, a sepa- 
rate random process termed S, which is restricted to nonnegative values, is defined. 
The value of the S process is the modulus of the value of the S process. The rela- 
tions between these two processes are developed in appendix A. 

The probability density function of the Press process is determined by equation (3) 
and by the functions of the two component processes 

= ( 12 ) 

This relation was first derived by Bullen (ref. 4). The indicated function is a modified 
Bessel function of the second kind and zero order. (See ref. 22.) The corresponding 
characteristic function is obtained either from the component processes by using equa- 
tion (7) or from the Fourier transformation of the probability density function (eq. (12)). 

Cz(d) = + l) (13) 


By using the asymptotic expansion for the mo(iified Bessel function, the probability 
density fvmction of the Press process has basically an exponential form for large values 
of the argument 


p^Clzl » bA) - 


1 

\/27rAb|z| 



(14) 


The strongly non-Gaussian nature of the Press process is evident from the 
exponential -tjrpe dependence of the probability density function in its extremes. Another 
measure of the Gaussian character is the flatness factor, the ratio of the fourth moment 
to the square of the second moment. 


M4 = 


E[z^(tl) 

[Ep(tl)])^ 


(15) 


The moments are given by equations (10), The flatness factor has a value of nine for the 
Press process in contrast to the value of three for the Gaussian process. Thus, the first 
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order probability density function and the flatness factor both indicate the strongly non- 
Gaussian nature of the Press process. 


Vector -scalar . - The development is easily extended to a vector random process. 
However, the vector process can be defined in several ways, depending upon whether the 
R process or the S process or both are considered to be vector processes. The pres- 
ent development is guided by previous applications of the Press process in the aeronauti- 
cal literature, which consider the vector nature of the process to be due solely to the R 
process 

zj(ti) = rj(ti) s(ti) (16) 


Thus, the vector process Z is formed by the product of a vector process R, whose 
components correspond to those of the Z process, and a scalar process S. The distri- 
bution of the R process is specified to be jointly Gaussian with appropriate covariance 
matrix 


E[ri(tl) rj(ti) 



(17) 


The joint probability density and characteristic functions are developed from equations (3) 
and (7), with appropriate modifications for the vector nature of the R process. The 
probability density function of the vector-scalar Press process is 


Pz(z) = 


2 

TT 


- 11/2 


1-m 


{ 2 tt )^ det (H) 



(18) 


i=l j=l 

The joint characteristic function is derived from the joint characteristic function of the 
R process by using the vector form of equation (7) 


m m 


Cz(P) “ ^ ^ ^i^ij^j 

i=l j = l 


- 1/2 


The quantity m is the dimension of the vector R and Z processes. 
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(19) 
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The various moments of the vector Press process can be obtained from the 
moments of the component processes by using the product relation and the independence 
of the R and S processes. 

The joint probability density function of the Press process determines the proba- 
bility density function of various combinations of the vector components. For an exam- 
ple, the sum of two components of a vector process is 

zsum = zi + Z2 = (ri + r2)s = rgumS (20) 

The random process Rsum is Gaussian since it is the sum of two Gaussian processes. 
Thus, the random process Za nm is a Press process since it is the product of two inde- 
pendent Gaussian random processes. Therefore, the probability density of the sum is 
that of a scalar Press process with appropriate variance, determined from equation (20). 
This example shows that a linear combination of Press processes, all having the same 
scalar S process, is itself a Press process. 

The vector Press process described is defined with a scalar S process. In this 
case the resulting components of the Press process are never independent. This prop- 
erty follows either f-rom the probability density function, from the characteristic function, 
or directly from the product relation (eq. (16)). 

Vector -vector.- A vector Press process can be defined in another manner as the 
product of a vector R and a vector S process: 

Zj(tl) = rj(ti) Sj(ti) (21) 

In this case the probability density and characteristic functions are not those of equa- 
tions (18) and (19). Also, in this case the components of the vector process can be 
independent. 


Generalized Press Process 

Distributions of the S process, other than the Gaussian, have been considered 
both in the original development (ref. 1) and in several other reports. An extensive list 
of possible S component processes and the resulting product processes has been com- 
piled by Houbolt (ref. 23). In the present section the choice of the S process is gen- 
eralized by consideration of a family of processes related to the Pearson type in family. 
This idea was originally suggested by Bullen (ref. 5), by Houbolt (ref. 23), and by 
Dutton and Thompson (ref. 6). 
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The S process is specified to be a member of a family of random processes 
related to the Pearson type HI family 

]^(f (ip) <" = 

The Gaussian R component process is not changed. The resulting family of random 
processes is called a generalized Press process or simply a generalized process of 
index n. The family includes the special cases of the Gaussian (n = 1/2) and Rayleigh 
(n = 1) distributions of the S component process. Thus, the Press process discussed 
in the previous section is a generalized process of index one-half. 

The probabilistic structure of the generalized process is determined from the rela- 
tions for the R and S component process. The specific relations for the moments 
and for both the probability density and characteristic functions are developed in 
appendix B. 


SECOND AND fflGHER ORDER DISTRIBUTIONS 

The time evolution of the Press process is considered in the present section. The 
appropriate properties are defined by the higher order distributions: the joint probabil- 
ity of the random process considered at several time points. The second order distribu- 
tions are most important since these are required for consideration of the derivative and 
the exceedance expression of the random process. Also, the second order distributions 
are related to the covariance and power spectral density functions. Both the R and S 
processes are specified to be stationary and to have zero mean values. All indicated 
derivatives of the random processes and all indicated statistical moments are assumed 
to exist. 


Product Process 

The product process is developed from the product of two independent scalar ran- 
dom processes. 

z(t) = r(t) s(t) (23) 

The autocovariance function is the joint moment of the random process at two time 

points 


E[z(ti) z(t 2 )] = ’i^z(|t2 - til) = 


(24) 
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By using the product relation, the autocovariance of the product process is determined by 
the autocovariance functions of the R and S processes 




(25) 


Since the three random processes are stationary, their autocovariance functions depend 
only upon the difference of the values of the two time points considered. 

The derivative of the product process is directly related to the derivatives of the 
component processes 

z(t) = r(t) s(t) + r(t) s(t) (26) 

The auto covariance function of the derivative of the product process can be developed 
either from equation (25) or from equation (26) as 

^zir) = -^z(t) = ^r(T) ^s(-r) + ^r(r) ^s(^) - 2 ^;(t) ^s(^) (27) 

Thus, the autocovariance function of the derivative of the product process cannot be 
directly related to those of the derivatives of the two component processes alone; there 
is a mixed term not identified with the derivatives of either component process. This 
term vanishes in the special case of the variance of the derivative, 

E[z2(t)] = ^f^(O) = E[r2(t)] E[s2(t)] + E[r2(t)] E[s2(t)] (28) 


The auto covariance functions of higher order derivatives can be developed in a similar 
manner. 

The relationship between the product and component processes can be developed in 
terms of their power spectral density functions. The three processes are assumed to be 
mean square continuous; thus, they can be expressed in a Fourier or spectral expansion. 
The associated power spectral density function is related to the autocovariance function 
by the Fourier transformation relations 


^2;(t) = ^ e^^^4>2;(w) dco 




( 29 ) 


I& 
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Corresponding relations can be written for the R and S processes, with appropriate 
changes in the notation. 

The power spectral density function of the product process is determined by the 
corresponding functions of the component processes. The product relation between the 
autocovariance functions (eq. (25)) transforms into a convolution relation between the 
power spectral density functions 


^z(^^) = y - w) ^s(^) 


(30) 


The relations between the autocorrelation and power spectral density functions have been 
given previously in reference 24. 

The probability density fimctions of the product process are determined from the 
corresponding functions of the component processes. The second order probability den- 
sity functions are considered, that is, the relations for the processes considered at two 
time points. The joint second order probability density fvinction of the product process 
and one of its components is determined by the product relation (eq. (23)) and the inde- 
pendence of the two component processes as 

P2s(zi,Z2;si,S2) =-j^^Pr(r 1,1*2) Ps(Sl,S2) (31) 

The second order probability density function of the product process is obtained from 
the joint density function as 


Pz(zbZ2)=j j |^^Pr(i*bi'2) Ps(si.S2) dsi ds2 (32) 

The corresponding second order characteristic function is obtained from the probability 
density function by using the properties of the associated Fourier transformations 


^OO ^OO 

Cz(01,^2) = \ \ Cr(0isi,02S2) Ps(si,S2) ^^2 

_QO _or> 


(33) 


The preceding three relations are similar to the corresponding first order relations 
(eqs. (2), (3), and (7)). The formulation can be extended to higher order distributions. 

It is sometimes useful to interpret the formulation in terms of the associated con- 
ditional process. The appropriate relations are developed either from equation (31) or 
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from equation (33) by using the definition of the conditional probability density function 


Pc(zi,Z2|si,S2) = Pr(ri,r2) (34) 

Cc(^1,^2|si,S 2) = Cr(0iSi,02S2) (35) 


The preceding development considers the product process in general terms. The 
subsequent development is restricted to specific forms of the two component processes. 


Press Process 


The Press process is a special case of the product process; both components, R 
and S, are specified to be stationary Gaussian processes. The scalar form of the Press 
process is considered first. The two scalar component processes are completely speci- 
fied by their covariance functions and zero mean values. The second and higher order 
distributions of the Press process are completely determined by the covariance functions 
of the components. The covariance functions can be expressed in terms of the autocor- 
relation functions and the variances as 


Ejr(ti) r(t 2 )J = 4fi.(r) = a 2 pj.(r) 
E s(ti) s(t 2 )] = 4's(t) = b^ Pg(T) 


z(ti) z(t 2 )] = ^'ziCr) = A^b^ p^ir) 


Pz('^) = Pr('^) Ps("^) 


(36) 


(37) 


The relation between the autocorrelation functions follows from the relation between the 
autocovariance functions (eq. (25)). 

Similar relations can be written for the first derivative of the random processes by 
using equation (27). The autocorrelation functions of the derivatives are determined from 
the autocorrelation functions of the random processes as 


Pz(^) 


P’z(^) 

Pz(0) 


(38) 
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The second order Gaussian probability density and characteristic functions of the 
two component processes are 




aj.^(T) = 1 - Pi.2(t) 


(39) 


Cr(^l»^2) = exp 



+ 2 Pr(r) 0 i6>2 + ^2 



(40) 


Similar relations can be written for the S process with appropriate changes in notation 
as indicated by equation (36). 

The second order probability density and characteristic functions of the Press pro- 
cess are obtained from the product relations (eqs. (31) to (35)) and the functional relations 
for the R and S component processes. However, in the subsequent analysis it is sim- 
pler to work with the characteristic functions. The second order characteristic function 
of the Press process is obtained by using equation (33) 


Cz(^l,02) = + 2pj,pg^i^2 + ^2^ + 



(41) 


i8i = Ab0i 


(i = 1, 2) 


The characteristic function must be inverted in both transform variables to give the sec- 
ond order probability density function of the Press process. This inversion appears to 
be intractable. However, some of the important properties of the process can be 
obtained from the characteristic function. The second order characteristic function 
defines the corresponding moments of the Press process, although these are more easily 
obtained directly from the product relation. It is noted that the joint characteristic func- 
tion depends only upon the difference of the values of the two time points, which is a con- 
sequence of the stationarity of the process. 

The second order characteristic function can be used to determine several proper- 
ties of the derivative of the Press process. The appropriate relations are given by 
Moyal (ref. 25). The characteristic function of the derivative is obtained directly from 
the second order characteristic function 
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( 42 ) 


r -1/2 

C^{9a) = lim Cz(4^, ^) = [l - ^dV;(0)] [l - /3d2p”(0)] 


- 1/2 


^d = Ab0d 


The characteristic function of the derivative is separated into two factors. Each factor 
is the characteristic function of a Press process by equation (13) and is dependent upon 
the derivative of only one of the component processes. By equation (26) the derivative is 
the sum of two processes, each of which is the product of two independent Gaussian pro- 
cesses, Thus, the derivative process is the sum of two independent. Press processes, 
which is the statement of equations (42). 

The joint characteristic function of the Press process and its first derivative is 
used to examine the independence of the random process and its derivative and to deter- 
mine the exceedance expression. The joint characteristic function is determined from 
the second order function as 


Czz(6z,0d) = f^z^ + [l - )3d^ Pr(0)] [l - )3d^ P’^(O)] 


- 1 / 2 ' 




(43) 


/3d = Ab0d 


Since the joint characteristic function does not separate into a product of two separate 
functions of the individual Fourier transform variables, the random process and its first 
derivative are not independent, although they are uncorrelated since the process is 
stationary. 

The joint probability density function of the random process and its first derivative 
is obtained by inversion of the joint characteristic function in both variables. The joint 
density function is used to determine the expected rate of exceedances of a given value of 
the process, often referred to as the exceedance expression (ref. 26), 


pOO 

N(z) = \ z Pzz(z,z) dz 


(44) 


Evaluation of both the joint probability density function and the exceedance expression 
from the joint characteristic function appears to be intractable. 
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The higher order distributions of the scalar random process are developed in a 
similar manner. Since both the R and S processes are Gaussian, their higher order 
probability density and characteristic functions are determined by the matrices of the 
autocovariance functions 


E[r(ti) r(tj)] = a 2 p^(ti - tjf 

> 

= ^V,ij 


(45) 


The matrix of the autocovariance function of the S process is defined by similar nota- 
tion. The higher order characteristic functions of the Press process are determined by 
the same procedure as that used for the second order function (eqs. (41)). The appro- 
priate relations are expressed in the forms for the higher order distribution of the sca- 
lar processes: the basic characteristic function relation (eq. (33)), the characteristic 
function of the R process, and the probability density function of the S process. The 
higher order characteristic functions of the Press process are then obtained by integrat- 
ing all the elements of the S process in the relation corresponding to equation (33) 


Cz 




mij = 


m 

f=l 


(46) 


The scalar Press process is formally identical to a vector Press process at a sin- 
gle time point if the time evolution of the S component is suppressed. The higher order 
characteristic function of a scalar Press process (eqs. (46)) can be reduced to the first 
order characteristic function of a vector Press process (eq. (19)). This relationship is 
shown by giving a unit value to all elements of Pg and by interpreting the values of the 
scalar R process at given time points as the components of a vector R process at a 
single time point. This relationship can be extended to develop the higher order charac- 
teristic functions of a vector process. These functions are obtained from equations (46) 
by replacing the variances and correlation matrices with the appropriate covariance 
matrices of the vector R and of the scalar or vector S component processes. 
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Generalized Press Process 


The generalized Press process is formed by the product of the Gaussian R pro- 
cess with a S process which is a member of a family of processes related to the 
Pearson type in family. The first order probability density functions of the S process 
family are given by equation (22), The second order functions of the generalized process 
are determined directly from the second order functions of the two component processes 
and the appropriate product relations (eqs. (32) or (33)). However, the higher order func- 
tions are more easily obtained in a different manner by using the relation between the 
generalized processes of various indices, that is, the generalized process of index 2n 
Is the sum of n independent and identically distributed Press processes. This relation 
(eq. (B9)) determines the higher order functions of the generalized process from those of 
the Press process itself. This approach is discussed further in appendix B, 

In summary, the probabilistic functions of the Press and the generalized Press pro- 
cesses are developed by using the representation of the Press process as the product of 
two independent random processes. The probabilistic structure of the Press process is 
completely determined from the given properties of the component processes. The char- 
acteristic functions of arbitrary order can be developed for both the scalar and vector 
Press processes by introducing the appropriate notation in equations (46). 

APPLICATIONS OF PRESS PROCESS 

The mathematical properties of the Press process are developed in the previous 
sections. The present section examines the relationship between the mathematical prop- 
erties and the intuitive concepts of the Press process as originally developed and as 
applied both to atmospheric turbulence measurement problems and to associated aircraft 
response problems. The intent of the present section is to show how the intuitive con- 
cepts are related to and largely justified by the mathematical development. Two general 
topics are considered. The first topic is the development of the approximate form of the 
random process which is used for the analysis of the dynamic response of linear systems. 
This approximation is based upon the concept of the S process being slowly varying so 
that it affects system response only in a static manner. This concept is termed the quasi - 
steady approximation. The second topic is the transition properties of the Press process 
and their relation to the stationarity and the Gaussian character of data modeled by the 
Press process. 


The Quasi -Steady Approximation 

The quasi -steady approximation of the Press process is discussed in the present 
section. The quasi -steady form of the process is developed from the exact random pro- 
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cess. The validity of the quasi -steady approximation as a method for the analysis of 
dynamic system response to a class of non-Gaussian random processes is discussed. 

The time variation of the Press process requires consideration of the time varia- 
tions of both the R and the S processes. The dynamics of the exact Press process 
are influenced by the dynamics of both component processes. However, the historical 
development (refs. 1 and 3) introduced a different concept which is an approximation of 
the exact process dynamics. This concept is based on the turbulence patch idea; the tur- 
bulence field consists of patches where the random turbulence velocity component (Z) 
varies with negligible change of the random intensity or magnitude (related to S). This 
concept introduces the idea of two time scales. On a local time scale the Gaussian R 
process varies at a constant value of the S process. On a global time scale the S 
process varies, and introduces the non-Gaussian nature of the product process. This 
concept is called the quasi -steady approximation. The corresponding random process is 
called the quasi -steady Press process or simply the quasi -steady process. 

The development and the limitations of the quasi -steady approximation are shown 
explicitly by the consideration of the response of a linear dynamic system. The system 
response or output is the convolution of the impulsive response of the system and the sys- 
tem input (refs. 19 and 20). 


Zout(t) = f h(r) Zin(t - t) dr (47) 

The input Press process is the product of two independent component processes (eq. (23)) 
so that equation (47) can be written as 


Zout(t) = r h(r) rin(t - r) sin(t - r) dr (48) 

It is assumed that the impulsive response of the system is essentially restricted to some 
time period T. In other words if the impulsive response is set to zero after a time 
period T, the change in the system response calculated by equation (47) is negligible. 
With this restriction the integration in equation (48) can be limited to the system response 
time period. 


Zout(t) ® \ h(r) rin(t - t) Sin(t - t) dr (49) 
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At this point the quasi -steady concept is introduced; the S component process is essen- 
tially constant over the system response time period T. With this assumption the S 
component can be removed from the integration since it is independent of the integration 
variable. 


T 

Zout(t) ~ Sin(t) r h(r) rin(t - t) dr (50) 

^0 

The integration now involves only the R component process. The integral is identified 
as the response of the system to the R component process. 

T 

J'out(t) = ^ h(r) ri„(t - t) dr (51) 

Since the input R component process is Gaussian, then the output R component pro- 
cess also is since it is a linear combination of Gaussian random variables. The system 
response is the product of two independent component processes and combines equa- 
tions (50) and (51) 

Zout(^) ^ ^out(t) Sin(t) (52) 

Thus, the probabilistic structure of the system output is completely determined within the 
accuracy of the quasi-steady approximation. Further, the system output and input have 
the same probabilistic structure since they are the product of a Gaussian R process and 
the same S process. 

The development of the system response (eq. (52)) shows the structure and assump- 
tions of the quasi -steady process. The R component process is always specified to be 
Gaussian. There is no restriction on the distribution of the S component process since 
this component is not affected by the system response. It is noted that this pattern of the 
restrictions on the two-component process is precisely the pattern used in the original 
development. (See refs. 1 and 3.) The time variation of the S component process is 
assumed to be much slower than all significant aspects of the system response. In appli- 
cations to atmospheric turbulence problems, the R component process is also assumed 
to vary on a much smaller time scale than the S component. This is the ^Tocally 
Gaussian’^ condition. The development of equation (52) shows that this assumption is 
actually not necessary. However, this assumption is followed in the subsequent develop- 
ment, since the case of a rapidly varying R component is of primary importance in 
dynamic response analysis. 


21 



The quasi -steady approximation presents a method of analysis of the response of 
dynamic systems to one class of strongly non-Gaussian processes. In applications the 
input R component process is defined to have unit variance; that is, 

Ejr?^(t)j = 1 (53) 

This definition establishes an arbitrary scale factor between the two component processes. 
The system response is determined by the response R component process. For linear 
systems, the response is completely specified in a probabilistic sense by the covariance 
function of the Gaussian R component and the variance of the S component process. 

The probabilistic relations for the quasi -steady process are obtained from those 
for the exact Press process. Either of two analytical procedures can be used. In the 
first procedure the relations for the quasi-steady process are obtained by developing the 
dynamics of the Gaussian R process, and then introducing the quasi -steady S process 
by means of consistency relation. This method is the standard approach used in the 
development and application of the Press process in the aeronautical literature. This 
procedure is the same as that for developing the first order relations of a vector Press 
process formed by the product of a vector R and a scalar S component (eqs. (18) 
and (19)). In the second procedure the relations for the quasi -steady process are obtained 
directly from the corresponding relations for the exact process by setting the local vari- 
ation of the S process to zero. The second method explicitly shows the development of 
the quasi-steady approximation. 

An example of the quasi -steady approximation is the development of the derivative 
of the product process (eq. (26)) 

z(t) = r(t) s(t) + r(t) s(t) 

The relative contributions of the two component processes to the variance of the deriva- 
tive are expressed in terms of the derivatives of the two autocorrelation functions. 

E[z2(t)] = A2b2[-p”(0) - p’^(O)] (54) 

Under the quasi -steady approximation, the time variation of the S process is negligible 
in comparison with that of the R process. The autocorrelation function of the S com- 
ponent is constant on a local time scale. The variance of the derivative of the quasi- 
steady process is 

E[x2(t)] = -A2b2p’^(0) (55) 
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(The notation of the Z and X processes is used to distinguish the exact and quasi- 
steady processes, respectively.) 

Similarly, the characteristic function of the derivative of the quasi -steady process 
is developed from that of the exact Press process given by equations (42) 


Cz(9) = [l 

= [l 


- 


Pil(O)] 

P’:^(0)] 


- 1/2 


1 - 


- 1/2 



(56) 


Thus, the derivative of the exact process is the sum of two independent Press processes. 
The derivative of the quasi-steady process is itself a Press process. This relation is 
equivalent to obtaining the derivative of the quasi -steady process from the product rela- 
tion by omitting the derivative of the S component 

x(t) = r(t) s(t) (57) 

One important application of the quasi -steady concept is the development of the 
associated exceedancd expression, which is the expected rate of crossings of a given level 
of a random process. The exceedance expression is obtained from the joint probability 
density of the random process and its first derivative (ref. 26) 


N(x) = y X Pjjjj(x,x) dx 


(58) 


An alternate approach is to use the exceedance expression of the R process which is 
Gaussian. The exceedance expression of the total process is obtained by using the con- 
sistency relation to eliminate the S process 


N(x) = \ N(x|s) ps(s) ds (59) 

The exceedance expression for the quasi -steady Press process is determined by either 
of these approaches. Historically, the development of the exceedance expression has 
followed equation (59) 


N(x) = No exp 


(60) 
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The exceedance expression for the quasi -steady generalized Press process of index n 
is obtained in a similar manner since only the distribution of the S process is changed 


N(x;n) 


Nq /|x|V 

gn-l r(n)W 



(61) 


The quasi -steady process is an approximate form of the exact Press process. The 
time variation of the S component process is assumed to be much slower than both the 
R component process and all significant aspects of the dynamic response of a system. 
The quasi -steady assumption is an approximate method for the analysis of the response 
of linear dynamic system to a class of strongly non-Gaussian random processes. 

Transition Properties 

The transition properties of the Press process are examined in the present section. 
By using the transition properties, the relation between the properties of the Press pro- 
cess and the properties of stationarity and Gaussianness is also examined. 

Measured exceedance data for atmospheric turbulence velocity components show a 
dual nature, that is, both Gaussian and exponential -type distributions. The Gaussian dis- 
tribution appears in some of the direct measurements of atmospheric turbulence, gener- 
ally those related to short periods of turbulence. The exponential -type distribution also 
appears in some of the direct turbulence measurements and is particularly valid for 
airplane VGH data measured over a long period of time. (See refs. 12 and 27.) A dis- 
cussion of the existence of both types of data is given in reference 28. 

The original development of the Press process (ref. 1) attempted to account for this 
dual nature by introducing the turbulence patch concept. Thus, data measured within one 
patch are Gaussian; data measured over a large number of patches have an exponential- 
type distribution. The transition between the two distributions is interpreted as an indi- 
cation of nonstationarity. 

The questions of the stationarity and the Gaussian character of the Press process 
are formally answered by the functional form of the probability density functions; the 
Press process is stationary and non-Gaussian. The formal definition of these properties 
does not answer the questions regarding experimental data. However, these questions 
can be restated in a more specific form. First, how can the Press process, which is 
stationary, account for data which appear to be nonstationary? Second, how can the 
Press process, which is non-Gaussian, account for data which appear to be Gaussian? 

The answers to both of these questions are contained in the concept of the product of two 
independent random processes, which evolve on radically different time scales. 
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The transition effects of a random process are described by a transition probability 
density fimction, the probability density fvinction of a random process conditional on its 
value at an earlier time point. Any measurement of a random function implicitly involves 
consideration of a transition probability density function which determines the transition 
of the process from a given initial state to a random final state. In other words, a set 
of measured data, started at some point in time, must forget its initial value to some 
extent. With the usual type of random process the transition effects are eliminated by 
measuring data over a time period which is many times greater than the time scale of 
the random process. With the Press process this operation is complicated by the pres- 
ence of two independent random processes that have radically different time scales. 

Thus, the transition process occurs in two steps and is indicated schematically by the 
following notation: 

p[z(t 2 )|r(ti),s(ti)] 

- p[z(t2)|s(ti)] 


=(*2) 


(ti < t2) 


The first transition involves the loss of the initial value of the R process. For the 
case of atmospheric turbulence, the time scale of this process is usually on the order of 
a few seconds, which is usually much less than the measurement period. Thus, this 
transition is of secondary importance in most applications; it is omitted in the subsequent 
discussion. The second transition, involving the loss of the initial value of the S pro- 
cess, may involve a time scale on the order of the measurement period. This transition 
may strongly influence the properties of the measured data. 

The important measured quantity is the conditional or transitional probability den- 
sity function of the random process at time t2, given the value of the S process at an 
earlier time ti. This function is determined from the consistency relation and the defi- 
nition of conditional probability as 


p(z2|si) 


P(z2|si) 



P(z2|si»S2) P(s2|si) ds2 




y P(Z2|S2) P(s2|si) dS2 

-oo 


(62) 
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The second relation follows from a special functional property of the conditional 
probability density function 

P(Z2|®1’S2) =P(Z2|®2) (63) 

This relationship is proven by using the properties of the associated characteristic 
function 


Cp2|si,S2) = Cc(0i = 0,e2|si,S2) 

> 

C(02|si,S2) = Crpisj = O,02S2) 


(64) 


The first equality is the consistency relation for the characteristic function which obtains 
the marginal fvinction from the second order characteristic function. The second equality, 
which is a functional property (eq. (35)) of the second order conditional characteristic 
function, shows that the indicated characteristic function is independent of s(tj). 

The relation (eq. (63)) does not imply an independence property between the product 
process Z and the S component process. 

The two probability density functions under the integral sign in equations (62) are 
both Gaussian for the Press process. The Indicated integration cannot be performed 
directly but requires a series expansion and 

oo 

£=0 


Integration term by term. 





0!s b 


Ois = \fi 


Ps" 


^ 0 
^ Ab 


The convergence of the series follows from comparison with the multiplication theorem 
of Bessel functions (ref. 22), The series converges for all finite values of the param- 
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eter Accordingly, the series converges for all values of the modulus of Pg(T) less 
than one. 

The transitional probability density function (eq. (65)) defines the transition of the 
Press process from the initial random state to the final random state. The initial and 
final states are determined by considering the limiting values of Ps{t) as the time dif- 
ference becomes zero and infinity, respectively. 




llm p(z2|Si) . ^ Ko 




' 2 | 


Ab 


( 66 ) 


(67) 


Equation (66) follows from the integral expression (eqs, (62)). The development requires 
consideration of the Dirac delta functional form of the transition probability density func- 
tion of the S process in the limit of small time differences. Equation (66) does not fol- 
low from the series expression (eq. (6 5)), since that series is not uniformly convergent. 
Equation (67) follows either from the integral expression or from the series expression. 

The transition probability density function is Gaussian only in the limit of zero 
time difference; the Press process itself is non-Gaussian. Also, the transition probabil- 
ity density function depends only upon the time difference. This property is a direct con- 
sequence of the stationarity of the Press process, although the transition itself is 
nonstationary. 

The transition properties of the Press process can be developed in terms of an 
associated transition random process. This process is defined by the transition proba- 
bility density function (eqs. (62)). Comparison of equations (62) and (5) shows that the 
transitional process is the product of two independent random processes, the fully devel- 
oped R process and the transitonal S process 

ztr(t) = r(t) str(t) (68) 

The transitional S component process is defined by the probability density function of 
s(t2) conditional on the value of s(ti). It is Gaussian with known mean and variance, 
both dependent on Pq(t). 
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The transition process generates an associated set of transition moments 


P(^2|si) dz2 


( 69 ) 


These moments can be calculated directly from the product relation (eq. (68)) 
E(z,y = E(r«) E(sy 


^oo 

E(str) = j S2^ p(s2jsi) ds2 




(70) 


The transition moments are defined by ensemble averages only, since the transition pro- 
cess is nonstationary. The relation between the ensemble averages and the associated 
time averages is discussed subsequently. 

The moments show the transition of the Press process in a concise form. The sec- 
ond and fourth moments are of primary interest: 


= A2[b2 + pg2(r) (sj2 . 

e(z^j.) = 3b^ + 6 Pg2(T) ^si2 _ b2j + p ^{ t ) ^sj'^ - 6sj^b2 + 3b^) 


(71) 


The moments show the transition through the autocorrelation function of the S process. 
The transition variance shows the development of the variance of the process from its 
given initial value to its final value which corresponds to the fully developed Press 
process 


E 



a2sj2 

A2b2 


(t - 0) 
(t - oo) 


(72) 


The transition of the probabilistic structure is also shown by the corresponding flatness 
factor 


E 


M4,tr('^) = 7 


[4r) 


K4) 


-l2 


(73) 
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The transition flatness factor is plotted in figure 2 as a function of the autocorrelation 
function of the S process. The transition flatness factor varies from the value of 
three for the initial Gaussian distribution to the value of nine for the fully developed 
Press distribution. The development of the transition process is strongly influenced by 
the transition of the variance of the S process, from its given initial value to the vari- 
ance of the fully developed process. The transition of the Press process is slower in 
terms of Pg for higher initial values of the S process. 

The development of transition moments can be extended to the case of the general- 
ized Press process. The moments are developed from the corresponding transition pro- 
cess, which is generated in the same manner as the generalized process itself. Accord- 
ingly, the transition process is the sum of mutually independent and identically distributed 
transition Press processes, defined by a relation similar to equation (B9). 



2n 

ytr(t;n) = ^ 

j=l 


(74) 


The transition moments of the generalized process can be expressed in terms of the 
moments of the Press process itself, which is a generalized process of index n equal 
to one -half 


E(y?r) = 2n e(z2J 


E(y^J = 2n e(zJJ + 6n(2n - 1) e(z2^) 


(75) 


In a similar manner the transition flatness factor of the generalized process is related 
to that of the Press process 


M4,tr(^;n) = ^[M4,tr(^;“ = 1/2) - 3] + 3 


(76) 


The transition flatness factors for the generalized processes are plotted in figure 3 as a 
function of the autocorrelation function of the S process. The transition flatness factor 
varies from the value of three for the initial Gaussian distribution to the value for the 
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fully developed process given by equation (B6). The transition factors also show the 
trend toward a Gaussian distribution as the index of the generalized process is increased. 
Only the special case where the initial value is equal to the standard deviation of the S 
process is considered; the additional effect of the transition of the variance (eqs. (75)) is 
not shown. 

The transition properties of the Press process can also be examined by means of 
the exceedance expression, which is a commonly measured quantity in atmospheric tur- 
bulence applications. The quasi -steady approximation is used to simplify the develop- 
ment. The exceedance expression for the transition process is developed in the same 
manner as the expression for the fully developed process (eq. (59)) 

N(x2|si)=J N(x 2|S2) p(s2|si) ds2 (77) 


The transition exceedance expression is the expected rate of exceedances of a specified 
level at time t 2 given the value of the S process at an earlier time tx. The proce- 
dure for the integration of equation (77) is similar to that for equations (62) 



The notation follows that of equation (65). 

An example of the transition exceedance expression is shown in figure 4. The 
exceedance expression shows the development of the process from the initial Gaussian 
distribution through intermediate states to the final ejqjonential form of the fully developed 
Press process. The expected rate of exceedances is strongly affected by the transition, 
particularly for large values of the random process. The transition exceedance expres- 
sion can have a predominantly exponential form when pg is still fairly large. The 
example of figure 4 is the special case where the given initial value of the S process is 
equal to its standard deviation. Thus, the example shows the transition of the probabilis- 
tic structure but does not show the additional effect of the transition of the variance. 

Figure 5 is an example of the transition exceedance expression which shows the effects 
of the transition of both the probabilistic structure and the variance. In this case the 
development of the variance increases the transition effects upon the exceedance expres- 
sion, particularly for large values of the random process. 

In the preceding development the statistical moments of the transition process are 
defined by ensemble averages. Since the transition process is nonstationary, the ensem- 
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ble averages are generally not equal to the time averages of a single record of the pro- 
cess. For a stationary and ergodic process, the two types of averages can be equated 
with reasonable accuracy if the period of time measurement is sufficiently large. This 
relation can possibly be applied to the transition process since that process is stationary 
on a local time scale. Thus, the ensemble averages of the transition process can be 
replaced by time averages if the variation of the S process is negligible both with 
respect to the development of the R process and with respect to the required time mea- 
surement period of the R process. Both of these conditions are contained in the con- 
cept of the quasi -steady process. 

The measurement of experimental data over a sufficiently long period of time is 
influenced by the development of the S process. The transition exceedance expression 
developed previously is essentially an instantaneous quantity, although the expression 
includes the full development of the R process. The time-averaged exceedance expres- 
sion is related to the instantaneous quantity defined previously. 


<^N(x " T It ^ ^ ^1 


dt 


(79) 


The time-averaged expression depends on the length of the measurement period. The 
calculation of the averaged exceedance function from the instantaneous function (eq. (79)) 
requires knowledge of the time dependence of the transition process, which is completely 
determined by the autocorrelation function of the S process. 

Transition properties are generated by the evolution of a random process. Transi- 
tion properties are possessed by all random processes, excluding the case of a totally 
static process. For the usual random process these properties represent the transition 
from a discrete distribution at a given initial value to the continuous distribution of the 
fully developed process. In the Press process the development occurs in two stages 
because of the different time scales of the two component processes. The first stage rep- 
resents the development of the R process with negligible development of the S pro- 
cess. The second stage, which is considered in the previous discussion, represents the 
development of the S process. During this stage, the probability density function devel- 
ops between two continuous distributions rather than the usual transition from a discrete 
to a continuous distribution. 

The Press process has a formal simularity to the anal 5 d;ical model of evolutionary 
spectra used for the analysis of nonstationary processes. The concepts of evolutionary 
spectra and oscillatory processes have been developed by Priestley (ref. 29). These 
concepts have been applied to the problem of aircraft response to atmospheric turbulence 
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by Howell and Lin (ref. 30), and Verdon and Steiner (ref, 31), The relationship between 
the Press and the oscillatory processes is discussed in appendix C, 

In summary, the properties of the Press process answer several questions regard- 
ing the stationarity and the Gaussian character of measured atmospheric turbulence data. 
The answers to these questions are contained in the concept of the product of two indepen- 
dent random processes which evolve on radically different time scales. The transition 
properties of the Press process generate a comprehensive model of the nonstationary 
aspects of atmospheric turbulence. The transition properties of the Press process 
account for both the Gaussian and the exponential forms of measured data and explicitly 
account for the transition between these two forms. 

APPLICATIONS TO ATMOSPHERIC TURBULENCE PROBLEMS 

The application of the Press and the generalized Press processes to atmospheric 
turbulence problems is considered in this section. Two basic subjects are considered: 
first, the application of the processes to measured turbulence data and second, the effect 
of the distribution of the S component process upon calculated exceedances. A specific 
example of an airplane mission analysis is developed. The relation between the distribu- 
tion of the S process and the rate of the exceedance of the extreme loads associated 
with the determination of structural strength is examined. 

Measured Atmospheric Turbulence Data 

The application of the Press process to atmospheric turbulence measurements 
requires the introduction of the concept of several types of turbulence. This concept was 
introduced in the original development (ref. 1) and has extensive experimental justifica- 
tion. The S component process is replaced by a conditional process, which is condi- 
tional on the type of turbulence. The probability density function of the modified S pro- 
cess is obtained by the consistency relation 


Ps(s) =^Pi P(s|i) 
i 


(80) 


The parameters Pj are the probabilities of occurrence of the ith type of turbulence and 
are generally referred to as the probability parameters. Each of the conditional proba- 
bility density functions (eq. (80)) has a separate standard deviation, whose associated 
parameters bi are generally referred to as the intensity parameters. 

Most developments consider two types of turbulence, that is, two terms in the sum- 
mation in equation (80). These are generally referred to as nonstorm (i = 1) and storm 
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(i = 2) turbulence, althoi;^h the specific identification is questionable. There have been 
suggestions for a third term to represent extreme turbulence. There is an additional 
term, omitted in the literature, which represents the probability of no turbulence. This 
term does not contribute to any of the measured quantities except for the zero values of 
the random process. However, this term is necessary to insure that the total probabil- 
ity is equal to one. 

The modified Press process is formed by the product of the Gaussian R process 
and the modified S process. For the modified Press process each conditional distri- 
bution (eq. (80)) is Gaussian, The resulting S component process is not Gaussian if 
there is more than one type of turbulence, although it is usually identified as Gaussian 
in the aeronautical literature. The probability density fimction and the exceedance 
expression for the modified Press process are obtained from the relations for the com- 
ponent processes and equations (3) and (59) 



The mean value of the random process is included in the preceding expressions. The 
exceedance expression (eq. (82)) has been the basis of extensive experimental study of 
atmospheric turbulence, for example, references 11 to 16. A review of the experimental 
determination of the associated atmospheric parameters is given in reference 12. 

The procedure can be extended to the generalized Press process by using the rela- 
tions for the generalized process development in appendix B. 



The Press process has an inherent dependence property which affects different 
components both of a vector process and of quantities formed from higher order distri- 
butions. This property is caused by different quantities being generated with a common 
S component process. The dependence property is not changed either by the specific 
form of the S process or by the quasi -steady approximation. This property introduces 
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some subtleties in the properties of the Press process which are important in applications 
to atmospheric turbulence. Two examples are considered. 

The Press process and its first derivative are not independent. They are uncorre- 
lated since the process is stationary. However, the independence property is absent 
since the two quantities have the same S component process. One consequence of this 
property is that the exceedance expression and the probability density have different func- 
tional forms (eqs. (81) and (82)). This is consistent with the development of the exceed- 
ance expression from the joint probability density function of the process and its first 
derivative (ref. 26). Consequently, it is necessary to consider the exceedance and the 
probability density functions separately since they cannot be obtained from each other. 

The two functions can be related in the limit of large values of the argument of the mod- 
ified Bessel functions (ref. 22) as follows: 


N(x;n) 

Px(x;n) 
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(84) 


Examination of some measured atmospheric turbulence data (ref. 32) which were analyzed 
by both types of fimctions shows a qualitative agreement with equations (84). 

Another example of the dependence property of the Press process is the associated 
vector process. In atmospheric turbulence applications the vector process is assumed 
to be formed by the product of a vector R process and a scalar S process, that is, the 
vector-scalar form of the process. The components of the Gaussian R process are 
independent if they are uncorrelated, as is assumed for isotropic turbulence. The com- 
ponents of the resulting Press process are also uncorrelated. However, they are not 
independent. This property does not affect the joint moments of second order between 
different components of the turbulence, but it does affect the joint moments of higher 
order. If the vector -vector form of the Press process is used, then the components of 
the resulting process can be independent. 

The exceedance ratio and the probability distribution function have been used as 
measures of the rate or probability of failure of a system. For linear systems these 
measures can be related to a single intensity factor because of the specific functional 
form of the exceedance ratio (eqs. (82) and (83)) and the probability distribution function 
(from eq. (81)). These are functions solely of the ratio of the incremental response level 
and the standard deviation of the response R process for a given set of atmospheric 
parameters 
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Thus, this turbulence or gust intensity parameter is directly related to a value of the 
probability of exceedance for a linear system. Also, the gust intensity parameter and 
the expected number of zero crossings are directly related to a value of the expected 
number of exceedances of a particular load level. It has been the basis of the develop- 
ment of structural strength criteria for aircraft response to atmospheric turbulence 
(refs. 14 to 16). 

Effect of Intensity Distribution Upon Calculated Exceedances 

The effects of the distribution of the intensity (S) component of the product process 
upon the calculated exceedances of aircraft response and loads are examined in this 
section. The cases of a Rayleigh and of one other distribution of the intensity process 
are considered. The development is divided into two parts. First, methods of deter- 
mining the atmospheric parameters for different intensity distributions are developed. 
Second, the effects of the distribution of the intensity process upon calculated exceedances 
are examined. The exceedances are examined by both the mission analysis and the 
design envelope techniques (ref. 14) which are used to define the required limit. load level. 
The mission analysis technique is applied to the Boeing 720B for which a complete set of 
mission data is available in reference 15. 

The generalized Press process is used to consider different distributions for the 
intensity process. Two specific cases are considered. The first is a Rayleigh distribu- 
tion, which is a generalized process of index one. The Rayleigh distribution of the inten- 
sity process has been suggested by some examinations of atmospheric turbulence data 
(for example, refs. 11 and 33). The second is a higher order distribution, of index four, 
which shows the effects of a more radical change in the distribution of the intensity pro- 
cess. Both cases are compared with the Press process which has a Gaussian distribu- 
tion of intensity and which is a generalized process of index one -half. Only the fully 
developed form of the random processes is considered since the analysis is concerned 
with exceedances measured over a long period of time. The present section uses the 
terminology of the aeronautical literature; the S component is called the intensity pro- 
cess, whose classification is based upon one type of turbulence only. 

The use of the generalized process requires the determination of the parameters of 
the intensity process (the P’s and b’s) for different indices of the process. Different sets 
of the parameters are needed since different indices of the generalized process have dif- 
ferent exceedance expressions (eq. (83)). The historical approach is to determine the 
parameters so that the calculated exceedances match the measured exceedance data. A 
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slightly different approach is used in the present study. A method is developed which 
modifies the available values of the parameters which are defined for the Press process. 
This approach is adequate to determine the general effects of the intensity distribution. 
This ^proach has two advantages. First, it avoids a reexamination of the basic experi- 
mental data. Second, it shows the effect of the intensity distribution on a comparative 
basis and avoids some open questions on the values of the original parameters for the 
Press process. Thus, any revision of the original parameters, such as those sv^gested 
in references 12 and 34, will have a minor effect on the present comparison. 

The present approach for determination of the atmospheric parameters is based 
upon the idea of matching measured exceedance data. However, the matching of exceed- 
ance expressions does not imply the matching of other quantities as the index of the gen- 
eralized process is changed. One study (ref. 16) shows the impossibility of matching 
both the calculated exceedances and the variance for different indices of the generalized 
process. The present approach is based upon matching the exceedance expressions since 
these represent the basic experimental data used to determine the original parameters. 

The parameters for the generalized process are determined by the following pro- 
cedure. The original parameters of the Press process are scaled by two factors. The 
same factors are used for all altitudes. The original set of parameters is that of refer- 
ences 14 and 15. The same set of parameters is used in aircraft strength criteria 
(ref. 18 and (with minor revision) ref. 17). The scaling factors are determined by 
matching the two exceedance curves in the range of three to four times the standard devi- 
ation of the original Press process. As a check, the new parameters are then used in a 
mission analysis of the Boeing 720B load factor (vertical acceleration of the center of 
gravity). The calculated exceedances are compared with those for the Press process to 
see whether they match in the range of 0.6g to 0.8g, which is the primary range of the 
measured data (refs. 12 and 27). The procedure considers only the storm (i = 2) terms. 
The same factors are applied to the nonstorm (i = 1) terms, although these can usually be 
ignored at the limit load levels of interest. 

Thus, the calculated exceedances of the generalized process must match those of 
the Press process if the parameters are properly defined. This does not say, however, 
that the two sets of calculated exceedances are identical. Differences can occur for two 
basic reasons. First, differences can occur at different response levels such as those 
associated with limit load levels which may involve an extrapolation from the range of 
the measured data. Second, differences can occur for quantities other than the load fac- 
tor which is the basis of the measured data. 

One approach to the development of structural criteria is the mission analysis 
technique. This approach is based upon the exceedances in average aircraft us£^e. The 
aircraft usage is modeled by a set of representative flight segments. The total exceed- 
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ances are the sum of the exceedances calculated for each segment, weighted by usage 
factors. A discussion of this technique is given in reference 14. 

The mission analysis technique is applied to the Boeing 720B. A mission analysis 
of the 720B, using the Press process, was part of the FAA study to develop the structural 
criteria. (See refs. 14 and 15.) The required response and mission data are given in 
reference 15. Three response quantities are examined. First, the load factor is exam- 
ined primarily to check the adjustment of the atmospheric parameters. The second quan- 
tity is the critical load for the vertical (symmetric) analysis, wing bending moment at 
wing eta station 0.33. The third quantity is the critical load for the lateral analysis, 
vertical -tail bending moment at elastic axis station 158. For the last quantity the yaw 
damper is assumed to be fully operative by following the procedure of reference 15. In 
the development of structural criteria based on the mission analysis technique, an exceed- 
ance level of 2.0 x 10”® exceedances per average hour is often used to specify the required 
level of limit load (refs. 14 and 17). 

The first case considered is a generalized process of index one, which has a 
Rayleigh distribution of the intensity process. The atmospheric parameters are adjusted 
by the procedure indicated previously 

Pi(n = 1) = 0.59Pi(n = 1/2) 


bi(n = 1) = 0.89bi(n = 1/2) 


E[x2(t;n = 1) = 0.935E x2(t;n = 1/2) 


The adjustment factors for the probability and intensity parameters and the resulting fac- 
tor for the variance, using equations (B2), are indicated. The largest change is in the 
probability parameter. The changes in the intensity parameter and in the variance of the 
total process are smaller. 

The results of the mission analysis for the generalized process of index one are 
shown in figures 6 to 8. The results are the exceedances of load level per hour of aver- 
age airplane usage plotted against the load level. The exceedances for the load factor 
(fig. 6) show close agreement for the two processes. This agreement indicates that the 
atmospheric parameters have been adjusted in a reasonable manner. The exceedances 
for the two loads (figs. 7 and 8) show essentially the same comparison with small differ- 
ences at limit load levels. In all three cases the Press process gives a conservative 
specification of limit load; that is, it gives the highest values for limit load level. 
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The second case considered is a generalized process of index four. The adjustment 
factors and the resulting factor for the variance, using equations (B2), of the process are 
as follows: 


Pi(n = 4) = 0.25Pi(n = 1/2) 


bi(n = 4) = 0.60bi(n = 1/2) 

E[x2(t;n = 4)] = 0.72E[x2(t;n = 1/2)] 

All the quantities are changed significantly. The largest change is again in the probabil- 
ity parameter. 

The results of the mission analysis for the generalized process of index four are 
shown in figures 6 to 8, The exceedances of the load factor (fig. 6) show reasonable 
agreement in the range of 0,6g to 0.8g and indicate a reasonable determination of the 
atmospheric parameters. However, significant differences occur at higher response 
levels. For all three response quantities the differences at limit load levels are around 
10 percent. In all three cases the Press process gives a conservative specification of 
limit load. 

A second approach to the development of structural criteria is the design envelope 
technique. This criteria is based upon a single flight condition. The criteria is the 
intensity factor (eq. (85)) which implicitly specifies an associated exceedance ratio. Fig- 
ure 9 shows the value of the intensity factor as a function of altitude for specific values 
of the exceedance ratio (eq. (83)) and for the three generalized processes under consider- 
ation. This method follows the procedure of reference 14. For reference, a value of 
1.2 X 10"® is suggested for the criteria exceedance ratio in reference 14. The values of 
figure 9 indicate the effects of the distribution of the intensity process upon the intensity 
factor U(j. For the generalized process of index one, the differences from the Press 
process are insignificant. For the generalized process of index four, the differences are 
small. In both cases the Press process gives a conservative specification of limit loads. 
It is noted that there is a modification of this approach that uses the exceedances rather 
than the exceedance ratio to specify limit load levels (ref, 16). The present conclusions 
also apply to that approach since the only difference is the expected number of zero 
crossings of the R process which is independent of the intensity process. 

The preceding discussion shows that the calculated exceedances are largely inde- 
pendent of the distribution of the intensity process unless data are available over a wide 
rai^e of the response. Thus, it is generally difficult to determine the distribution of the 
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intensity process from measured exceedance data. A better approach is the use of the 
flatness factor, which is more sensitive to the distribution of the intensity process. If 
the measured data represent only one type of turbulence, the flatness factor is directly 
related to the index of the generalized process. For a fully developed process, the flat- 
ness factor has a simple form 




( 86 ) 


In summary, the effects of the distribution of the intensity process upon calculated 
exceedances are examined by using the generalized process. The exceedances are com- 
pared with those of the Press process which has a Gaussian intensity distribution. Sev- 
eral conclusions are made. First, it is necessary to redefine the atmospheric param- 
eters (P’s and b’s) as the intensity process is changed. The parameters can be defined 
so that the calculated exceedances reasonably, match the measured data. Second, the use 
of different distributions for the intensity process has a minor effect on the limit loads 
determined by the mission analysis technique. For a Rayleigh intensity distribution the 
effect is negligible. For higher index distributions the differences may not be negligible. 
However, the Press process, with a Gaussian intensity distribution, gives a conservative 
specification of limit load levels. Finally, it is generally difficult to determine the dis- 
tribution of the intensity process from measured exceedance data. A better method is 
the examination of flatness factors to determine the corresponding intensity distribution. 

CONCLUDING REMARKS 


The present report is an examination of the random process, termed the Press pro- 
cess, which is used to model atmospheric turbulence in aircraft response problems. The 
primary conclusion of the present study is that the basic concepts of the Press process, 
which were originally developed largely on an intuitive basis, can be derived from an 
explicit mathematical model. The model is based upon the product of two independent 
random processes, which have radically different time scales. With this model it is 
possible to describe explicitly the general properties of the Press process such as the 
locally Gaussian condition, the Gaussian and non-Gaussian forms of experimental data, 
and the possible nonstationarity of experimental exceedance data. 

The concept of the product process introduces two independent component processes. 
One component is the R process^ always specified as Gaussian, which defines the vari- 
ation of the product process at constant intensity or standard deviation. The other com- 
ponent is the S process which has a slower variation. The first, second, and higher 
order probability density and characteristic functions are obtained for the product process. 
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In developing these relations, both the usual Gaussian process and a general class of pro- 
cesses, related to the Pearson type III family, are considered for the S component 
process. 

The concept of a product process leads to an approximate procedure for the analysis 
of the response of linear djmamic systems to a class of non -Gaussian random processes. 
The procedure is based upon the concept that the time variation of the S component pro- 
cess is much slower than that of both the R component process and all significant 
aspects of the system d3mamics. Under this quasi-steady assumption the S3^em dynamics 
are influenced only by the R process, which is Gaussian. The variation of the S com- 
ponent process does not affect the system dynamics, but does introduce the non-Gaussian 
nature of the random process in a static manner. 

In the mathematical model the random process develops in two stages since the 
process is the product of two independent processes which develop on radically different 
time scales. The second stage of development is of primary interest, the development 
of the slower S component process with the R component process being fully devel- 
oped, The mathematical model explicitly shows the transition from the Gaussian form of 
the R component process to the non-Gaussian form of the fully developed Press process. 
This transition, which is nonstationary, is e?q)licitly described by the transition properties 
of the stationary Press process. 

The effects of the form of the distribution of the S component process upon cal- 
culated exceedances are examined. The cases of a Rayleigh and of a second distribution 
are compared with the Gaussian case. It is concluded that the calculated exceedances are 
not greatly influenced by the distribution of the S component if the atmospheric param- 
eters are properly defined, A specific numerical example of a mission analysis is com- 
puted for one airplane. The limit load levels determined by the criteria specification are 
compared with those computed from the Press process, with a Gaussian S component 
process. The changes due to the Rayleigh distribution for the S process are negligible. 
The use of the second process causes some differences at limit load levels. However, in 
the cases considered, comparison with the Press model shows that the differences are 
either negligible or that the Press model predicts somewhat higher occurrence of limit 
loads; that is, the Press model is conservative. 

The present study outlines the anals^tical properties of the Press process. Several 
applications of the formulation to the measured properties of atmospheric turbulence and 
several extensions of the analysis are suggested. 

Several properties of measured turbulence data which have not been examined on an 
explicit basis are identified. The primary areas of interest are the transition properties 
and the associated nonstationary process. The Press model predicts explicit relations 
for the transition properties of the exceedances, probability density function, moments. 
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flatness factor, etc. This prediction raises the basic question of the correspondence 
between the predicted properties and the measured properties of atmospheric turbulence: 
does the anal 5 rtical model adequately match the nonstationary properties of the measured 
data? An associated question is the specific form of the time variation of the S com- 
ponent process: what are the power spectral or covariance functions of the S process? 

The present study suggests two problems of a mathematical nature. First, the use 
of the transition properties of a stationary product process presents a method for the 
analysis of nonstationary random processes. Some of the properties of the associated 
nonstationary process are outlined in the present study. However, the relation between 
the present approach and other approaches to the analysis of nonstationary processes 
presents an interesting mathematical study. Second, the concepts of the product process 
and of the quasi -steady approximation present a method of analyzing the response of lin- 
ear dynamic systems to a class of non-Gaussian random processes. The justification 
for the approximation is discussed largely on a qualitative basis in the present study. 

This procedure should be thoroughly examined as a general analytical method, 

Lai^ley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
July 22, 1975 
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APPENDIX A 


DEVELOPMENT OF S PROCESS 


There are two random processes associated with the S component of the Press 
process. One is the S process, which is associated with the historical development 
of the Press model. The other is the S process itself, which is required for a com- 
plete mathematical formulation. The relations between these two processes are devel- 
oped in the present section. 

The difference between the two random processes is the range of their values: the 
S process is restricted to positive values whereas the S process has both positive and 
negative values. The values of the S process are the moduli of the values of the S 
process 

a = |s| (Al) 


The S process appears when the Gaussian R component process is related to 
a conditional process by using equations (4) and (11), 


P„(z|s) = —=i — 
\/27rAa 


ejq) 


-z‘ 


.2A^a2j 


(A2) 


The conditional process is Gaussian with standard deviation equal to Act which must be 
positive. Thus, the historical development, using the idea of a conditional process 
through the locally Gaussian condition, is based on the S process. Also, the analysis 
of turbulence data in terms of its standard deviation leads to the S rather than the S 
process. 

The properties of the S process are determined from the specified properties of 
the S process by equation (Al). The properties are obtained either from the probabil- 
ity density function of the S process or directly from the properties of the S process. 
For example, the moments of the S process are the absolute moments of the S pro- 
cess. The relations between the dynamic properties of the two processes are more com- 
plicated. The quantity of primary interest is either the autocovariance or the autocor- 
relation function of the E process 
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For the case of the Gaussian S component process, the relation between the autocorre- 
lation functions is (ref. 19) 

P(y{T) = PgiT) + |(sin r 

cos y = Pg(r) 

The dynamic properties of the two processes are determined from their autocorre- 
lation functions. The specific relations are obtained from equations (A4) either by eval- 
uating the derivatives of the autocorrelation functions or by expanding the functions in 
power series in the time difference variable. These procedures give the following 
relations: 



P'a(0) = 0 


p'iiO) = p”(0) 


p---(0+) = ^[-P’’(0)] 


3/2 


(A5) 


Thus, the variances of the first derivatives of the S and S process are equal. How- 
ever, the S process does not possess a second derivative since the corresponding var- 
iance has an infinite value. This is the primary reason the S process must be used in 
the formulation of the Press process: its dynamic properties can be specified arbitrarily. 
Thus, the quasi-steady approximation cannot be developed if the D process is used, 
since the second derivative of that process does not exist. However, once the quasi - 
steady approximation is established, the S process can be used in place of the S pro- 
cess since only the first order properties are used. 

In the aeronautical literature the S process is often referred to as the intensity” 
of the turbulence. This interpretation is based upon the standard deviation of the condi- 
tional process (eq. (A2)) being equal to Aa. Strictly speaking, the term is a misnomer 
since an intensity is related to the variance and not the standard deviation of a random 
process. However, the term intensity is commonly used to identify the S process. The 
S process is also identified as the intensity, if it is remembered that the standard devia- 
tion of the conditional process depends on the modulus of the S process. 

In summary, both the S and S processes are important aspects of the associated 
component of the Press process. The S process, with both positive and negative values, 
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must be used in the mathematical formulation. The S process, restricted to positive 
values, is required for the analysis of turbulence data. Also, the properties of the S 
process cannot be completely determined from turbulence data through the associated 
conditional process (eq. (A2)). The properties of the S process can be developed from 
the specified properties of the S process by using the transformation of equation (Al), 
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APPENDIX B 


DEVELOPMENT OF GENERALIZED PRESS PROCESS 

The generalized Press process is formed by the product of the Gaussian R pro- 
cess and the S process, which is specified as a member of a family related to the 
Pearson type in family of processes. The probability density functions of the S pro- 
cess family are given in equation (22). The moments and both the probability density and 
characteristic functions of the corresponding generalized processes are developed in this 
appendix. The relation between generalized processes of various indices is also 
develcped. 


First Order Distributions 

The variance of the scalar S process is determined from its probability density 
function (eq. (22)) 


E 




= 2n b^(n) = 2nb^ 


(Bl) 


The parameter b is not equal to the standard deviation of the S process in the general 
case. Also, this parameter can have different values for different generalized processes, 
that is, for different values of the index n, although the specific notation is dropped in 
the subsequent discussion. 


The moments of the scalar product process are determined by the moments of the 
R and S component processes and by the product relation (eq. (1)). (The generalized 
process is denoted as the Y process to distinguish it from the Press process itself.) 


y^(ti;n) 


= E 


-(tl) 


■(ti;n) 





in = odd) y (B2) 


(jf = even) 

J 
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The probability density and characteristic functions of the scalar generalized pro- 
cess follow from the explicit forms of the functions of the R and S processes and 
from application of equations (3) and (6) 



The formulation can be extended to vector processes without difficulty. The structure of 
the vector -scalar process is expressed most concisely by the characteristic function, 
which has a functional form similar to that of equation (19) 


( m m ' 

I I + 1 

1=1 i=i 


(B5) 


The generalized Press process has a strongly non -Gaussian character, as indicated 
by the exponential -type dependence of the probability density functions for extreme values 
of the process. Another indication is the associated flatness factor 


M4 



3(n + 1) 
n 


(B6) 


The values of the flatness factor are larger than the value of three for a Gaussian process. 
However, the flatness factor does approach the Gaussian value of three in the limit of 
large values of the index n. 

The index of the generalized process appears only as a power coefficient in the 
characteristic function (eq. (B5)). This property indicates a relation between the gener- 
alized processes of various indices, since the characteristic function of a sum of inde- 
pendent random variables is the product of the individual characteristic functions. The 
characteristic functions of the generalized process (eq. (B4)) and the Press process 
(eq. (13)) have a simple functional relation 
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Cy{6 ;n) = [Cy(0 ;n = 1/2)] = [C^ (0)] (B7) 

Thus, the sum of 2n independent and Identically distributed Press processes is a gen- 
eralized process of index equal to n 

2n 

y(ti;n) = 2 Zj(tl) (B8) 

j=l 

This relation presents an alternate method for developing the properties of the gener- 
alized process. For example, the relations for the moments (eqs. (Bl) and (B2)) can be 
develcped by this method. 

The development suggests a relationship to the central limit theorem. Interpreting 
the generalized process of index n as the sum of 2n independent and identically dis- 
tributed Press processes, the central limit theorem requires that the generalized process 
becomes Gaussian in the limit of large n. The validity of this statement can be shown 
by examination of thq characteristic function (eq. (B4)). This conclusion explains why 
the flatness factor (eq. (B6)) approaches the Gaussian value of three in the limit of 
large n. 


Second Order Distributions 

The development of second and higher order distributions of the generalized Press 
process requires specification of the higher order distributions of the S process. For 
the Press process, the S process is stationary and Gaussian; the higher order distri- 
butions are specified directly. For the generalized process it is simpler to use an alter- 
nate approach. 

For the first order distributions, the generalized Press process can be developed 
from the Press process itself by equation (B8). This relation is extended to the case of 
random processes considered at multiple time points 

2n 

y(t;n) = ^ Zj(t) (B9) 

i=i 

The random process Y is a generalized Press process of index n if the Zj pro- 
cesses are all Press processes (generalized processes of index one -half), which are 
independent and identically distributed. By using the relationship between the gener- 
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alized process and the Press process (eq. (B9)), the second order characteristic functions 
are developed 


Cy{0i,e2;n) = [Cz(0i,02)]^“ 


(BIO) 


This relation is similar to equation (B7) for the first order distribution. The second 
order characteristic function can be used to develop relations for the first derivative, the 
joint distribution of the process and its first derivative, and the exceedance expression. 
Thus, the second and higher order distributions for the generalized processes of various 
indices can be developed directly from the corresponding distributions of the Press 
process. 
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COMPARISON OF THE PRODUCT AND OSCILLATORY RANDOM PROCESSES 

The relation between the definition of the product process and that of the oscillatory 
process developed by Priestley (ref. 29) is examined in this appendix. It is shown that the 
spectral expansions of two random processes have a formal similarity. The discussion 
is restricted to the spectral expansions; the probabilistic structures of the random pro- 
cesses are not considered. 

The time variation of a random process is described by its spectral expansion. The 
product random process and its two component processes are stationary and are assumed 
to be mean square continuous. Accordingly, they can be expressed in a Fourier or spec- 
tral expansion (refs. 25 and 35) 





(Cl) 


The process {dR(cti)^ is the spectral process associated with the R process. The spec- 
tral process, having orthogonal increments, is the basis of the definition of the spectral 
distribution and power spectral density functions of the original process. 


E 



= dFp(ct)) = 


(C2) 


Similar relations can be written for the S and Z processes. By using the product 
relation, the power spectral density functions of the three processes can be related and 
will lead to equation (30). 

The Fourier expansion of the product process offers some insight into the time 
evolution of the process and into the relationship with the concept of a nonstationary pro- 
cess. By using the interpretation of the S process as slowly varying in time, the spec- 
tral expansion of the product process is 


z(t) = r(t) s(t) = j e^^* s(t) dR(o)) 


(C3) 


This relation follows directly from the product relation (eq. (23)) and the Fourier expan- 
sion of the R process (eq. (Cl)). 
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A similar relation is written for the Fourier expansion of a nonstationary process 
by following the concepts of evolutionary spectra and oscillatory processes developed by 
Priestley (ref. 29) 

z(t) = y ao(t) dZ(a>) (C4) 

-oo 


This relation is the spectral expansion of an oscillatory process in the special case of a 
uniformly modulated process. The modulation function ao(t) is a slowly varying 
deterministic function. 

The spectral expansions of the oscillatory process (eq. (C4)) and the product pro- 
cess (eq. (C3)) have a strong formal similarity. In both cases a spectral process is 
modulated by a time function which is considered to be slowly varying relative to the 
unmodulated process. With the oscillatory process the modulation function is a deter- 
ministic function; the resulting oscillatory random process is nonstationary. With the 
product process the modulation, function is a random process; the resulting product ran- 
dom process is stationary. The stationarity of the product process is determined by the 
functional properties of the probability density and characteristic functions and by the 
definition of stationarity. The nonstationary aspect of the product process requires an 
examination of the transition properties, which are considered in the main text. 
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Figure 1,- Atmospheric turbulence model with locally Gaussian regions 

and intensity variations. 
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Figure 3.- Transition flatness factor for generalized Press process 
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Figure 8.- Mission analysis for vertical -tail bending moment, generalized Press process. 
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